Survival of the black hole's Cauchy horizon under non-compact perturbations 
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We study numerically the evolution of spactime, and in particular of a spacetime singularity, inside 
a black hole under a class of perturbations of non-compact support. We use a very simplified toy 
model of a spherical charged black hole which is perturbed nonlinearly by a self-gravitating, spherical 
scalar field. The latter grows logarithmically with advanced time along an outgoing characteristic 
hypersurface. We find that for that class of perturbations a portion of the Cauchy horizon survives 
as a non-central, null singularity. 
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Introduction and summary: The geometrical and phys- 
ical properties of the Cauchy horizon singularity inside 
black holes have received much attention PJ . That singu- 
larity was shown to be null, non-central, and weak. The 
weak nature of the Cauchy horizon singularity has far- 
reaching implications. In particular, it leaves open the 
possibility that physical objects which fall into a black 
hole may traverse the Cauchy horizon singularity only 
mildly affected, and re-emerge in another universe. 

The evolution of spacetime geometry into a (weak) cur- 
vature singularity at the Cauchy horizon has been studied 
both numerically and analytically J|, ||[ ||, [|, 0, 0, ||, |j| . In 
all these studies the black hole was taken to be isolated, 
and the source of the perturbations was taken to be the 
perturbations which result from the evolution of nonvan- 
ishing multipole moments during the collapse. These per- 
turbations are inherent to any nonspherical gravitational 
collapse, and result from the backscattering of waves, 
which arc created during the collapse, off the curvature 
of spacetime |10 . These perturbations have a compact 



support at some initial time. 

It is interesting to ask whether the evolution of a null 
and weak singularity at the Cauchy horizon is just an 
artifact of the assumption of compactness. That is, 
will any dominating perturbation field which has non- 
compact support on the initial time slice lead to the full 
destruction of the null singularity, and its replacement 
by a spacelike one? In fact, the perturbations due to the 
collapse and the resulting tails can be thought of as a 
lower bound on the perturbation field. It is intersting to 
ask what happens to the Cauchy horizon if perturbations 
which are stronger than that lower bound are present. 

This question is interesting not just from the math- 
ematical viewpoint: indeed, a generic class of pertur- 
bations exists, where the perturbation field has non- 
compact support. These are the perturbations which 
arise from the capture of photons which originate from 
the relic cosmic background radiation (CBR). Even if 
removed from any conceivable astrophysical object, any 
black hole is still perturbed by the CBR. Unlike the 
perturbations due to the collapse, the CBR perturba- 
tions are non-compact. Because the perturbation field is 
greater than the lower bound set by the perturbations 
due to the tails, one may suspect that the evolution of 



spacetime, and in particular of the singularity, be domi- 
nated by the non-compact perturbations, rather than by 
the compact ones. 

If this is indeed the case, the non-compact perturba- 
tions threaten to change our notions of the causal struc- 
ture inside black holes. When the perturbations due to 
the tails are considered, it is found that the weakness 
of the Cauchy horizon singularity is a rather delicate 
issue: it depends on certain integrals being bounded. 
Specifically, assume that the field is due to the tails. 
Then, on the event horizon the scalar field <I> behaves 
like <& = (nv)~ n , where k is a constant, v is advanced 
time, and n is a positive integer which is related to 
the multipole moment of the perturbation field. De- 
noting schematically by TZ the fastest growing compo- 
nents of the Riemann-Christoffel tensor approaching the 
Cauchy horizon, the curvature at that limit behaves like 
TZ ~ t~ 2 [— ln(— r)] 2n+2 , where r is proper time along a 
timelike geodesic which is set equal to zero on the Cauchy 
horizon. The Cauchy horizon singularity is weak, if TZ is 
twice-integrable. (This last statement can be made pre- 
cise jo].) For positive value of n this is indeed the case. 
However, small changes in TZ, e.g., another factor of r~ e 
for any small and positive e, would change the picture en- 
tirely, as TZ would no longer be twice integrable. That is, 
the twice integrability of TZ is strongly dependent on the 
form of the field at the event horizon. The twice integra- 
bility of TZ (and consequently the weakness of the Cauchy 
horizon singularity) depends then on the assumption that 
the scalar field has a compact support (as this condition 
leads, through Price's analysis [^O), to the tail form for 
the field on the event horizon) . Dominating non-compact 
perturbations threaten to change TZ in a significant way, 
such that it would no longer be twice integrable. One can 
ask then the following question: Is the Cauchy horizon 
necessarily utterly destroyed and replaced by a spacelike 
singularity when perturbations with non-compact sup- 
port are present, or can it still survive (as a null, weakly 
singular hypersurface) also when perturbations of non- 
compact support are present? 

In this paper we shall answer the latter question in the 
affirmative. We show that a certain class of non-compact 
perturbations still preserves the null, non-central nature 
of the Cauchy horizon singularity, even though the evo- 
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lution of geometry and of the singularity is indeed domi- 
nated by the non-compact perturbations. We emphasize 
that any perturbation field with non-compact support 
(whose dynamics dominates deep inside the black hole) 
is appropriate, as it serves as a counter-example for the 
claim that only perturbations with compact support can 
evolve into a null and weak singularity. The survival of 
the Cauchy horizon as a null, non-central singularity can 
occur also when perturbations of non-compact support 
(of certain classes) are present. 

Model: In this paper we study the evolution of space- 
time curvature inside a black hole in the presence of 
perturbations which have non-compact support under a 
very simplified toy model. For simplicity, we take the 
black hole to be spherically symmetric, and to have a 
fixed electric charge Q. This is a useful toy model for a 
spinning black hole, because the unperturbed spacetimes, 
namely the Reissner-Nordstrom and Kerr spacetimes, re- 
spectively, have very similar causal structures, which lead 
to similar blue-sheet effects near their inner horizons. 
In fact, much of the understanding we currently have 
about black hole interiors have been obtained through 
the study of spherical charged models. (One important 
difference is that the null singularity inside a spherical 
charged black hole is monotonic, whereas the one inside 
a spinning black hole is oscillatory This difference 
is not crucial for our purposes here. Another difference 
is related to the question of the occurrence of a space- 
like singularity inside black holes. A spacelike singular- 
ity to the future of the Cauchy horizon singularity was 
found in spherically-symmetric, charged models. It has 
been argued that no corresponding spacelike singularity 
is likely to occur inside a rotating black hole [|l3|. Oth- 
ers have argued, that a spacelike singularity, possibly of 
the Belinskii-Khalatnikov-Lifshitz type, is a possible out- 
come. While this open question is extremely important, 
it is unrelated to the nature of the null singularity which 
precedes the spacelike one, if such a spacelike singularity 
exists.) 

We write the spherically-symmetric metric in double- 
null coordinates in the form 

ds 2 = -2e 2a ^' v) du dv + r 2 (u, v) dtf (I) 

where d£l 2 is the line element on the unit two-sphere. 
As the source term for the Einstein equations, we take 
the contributions of both the scalar field $ and the 
(sourceless) spherical electric field (see Q for details). 
The dynamical equations are the scalar field equation 
VpV M $ = and the Einstein equations, which reduce to 

*™« + -(r, u *„+r t ,$ u ) = (2) 
r 

r. u r, e 2 ° ( Q 2 \ 

and 

^ = 0. (4) 



These equations are supplemented by the two constraint 
equations 

r,uu - 2o-, u r, u +r($ u ) 2 = (5) 

r. vv - 2<7„r„ +r($„) 2 = 0. (6) 

Although similar, there is an important difference be- 
tween the numerical evolution of a code based on 
Eqs. (§)-(§) and a code which is based on the dynam- 
ical equations used in Ref. In the latter case, the 
wave equation for the g uv metric function becomes a free 
wave equation asymptotically close to the Cauchy hori- 
zon. (Notice that g uv vanishes exponentially in v near the 
Cauchy horizon, whereas r u ,r t ,4 tt and $ v decay like 
inverse powers of v, where v is proportional to advanced 
time - see below.) This implies that the numerical in- 
tegration becomes inaccurate near the Cauchy horizon 
because dynamically-important terms become negligible. 
When the field equations are written as Eqs. (§)-(§) this 
problem does not occur. 

Initial value problem: From the pure initial-value view- 
point, we need to specify three initial functions on each 
segment of the initial surface: r, a, and $. The constraint 
equations reduce this number: Eqs. (||) and (^) impose 
one constraint each on the initial data at u = Ui and 
v = Vi, respectively. The remaining two initial functions, 
however, represent only one physical degree of freedom: 
The other degree of freedom expresses nothing but the 
gauge freedom associated with the arbitrary coordinate 
transformation u — > u(u) , v — > v(v). In what follows we 
shall use a standard gauge, in which r is linear with v or 
u, correspondingly, on the two initial null segments. On 
the outgoing segment we take r. v — 1. (Notice, that this 
implies that this v is twice advanced time at late times.) 
On the ingoing segment, we take r u = const = r u o- 
(Notice that r U Q < 0.) The initial values of r are thus 
uniquely determined by the parameter ro = r(iii,Vi). We 
choose Ui — and Vi — ro, and thus we find: r v (v) = v, 
and r u (u) = tq + ur U Q. [Hereafter, we denote the ini- 
tial values of the three fields r, a, $ on the two segments 
of the characteristic hypersurface by r v (v),o~ v (v),<!> v (v) 
and r u (u), <r u (u), $ u (m), correspondingly.] Then, we can 
freely specify $ M (u) and $> v (v) (this choice represents a 
true physical degree of freedom) . The initial value of a is 
now determined from the constraint equations, namely 

o-u.u = r u ($u,u) 2 /(2r u0 ) , a v . v = r v ($ v . v ) 2 /2 , (7) 

together with the choice a(ui,Vi) = —(1/2) ln2. Thus, in 
the gauge we use, we need to specify two functions of one 
variable [<$> v (v) and and two parameters (ro and 

r u o) (in addition to the charge Q) for the characteristic 
initial value problem. 

Determination of characteristic data: We 
take the characteristic data to satisfy <& v (v) = 
(A/ yj —2r u o) \n(v/vi) along u = Uj. Here, A is a 
real constant, which is related to the amplitude of 
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the perturbation field. This choice for the scalar field 
is clearly of noncompact support, as the field grows 
logarithmically in advanced time. (Notice, that this 
implies that spacetime is not asymptotically flat.) In 
addition, we require that = on v = Uj, such 

that the field does not propagate outside the event 
horizon on the ingoing segment of the characteristic 
hypersurface. We also require that <E> is continuous at 

(tk,Vi). 

We note that is it unimportant what the field is along 
v = Vi and on u = m for Vi < v < Vf for any finite 
Vf. Any perturbation field with compact support leads 
to power-law tails at late times regardless of the spe- 
cific shape of the initial data. It is only the contribu- 
tions of the characteristic initial data from late advanced 
times which are important. That is, one can approxi- 
mate the characteristic hypersurface by dividing it into 
two parts: a compact part which is extended from (itj, Vi) 
to a point (ui,Vf) (with Vf > Vi), and a non-compact 
part which extends from (ui,Vf) forward, i.e., the points 
(iii,v > Vf). The specific form of the characteristic ini- 
tial data on (ui,v < Vf) is unimportant: it is only the 
contribution of the initial data at (ui,v > Vf) which is 
important. (Similarly, also the initial data along v — Vi is 
unimportant.) Consequently, we can determine arbitrary 
initial data at early times. 

The solution of the characteristic initial value prob- 
lem then is given by $ v (v) = (A/y/—2r u o) ln(v/vi), 
$„(«) = 0, a v (v) = -H2(v/v Q ) A2 /( 2r -o)]/2, and 
a u {u) = -(1/2) hi2. 

Numerical simulations: Our numerical code is a free 
evolution code in (1 + 1)-D in double-null coordinates 
with an adaptive mesh refinement (m). We tested the 
code and found that it is stable, and converges with sec- 
ond order. In the following we present results with the 
following choice of parameters, unless stated otherwise: 
-^initial = 1, Q = 0.95, tq = 5 (these uniquely determine 
the value of r u0 ), A — 0.3, and N = 10. Here, Mj n iti a i 
is the initial mass of the black hole, and N is defined 
as the number density of grid points on the character- 
istic hypersurface in both u and v directions. We find 
similar qualitative results also for other choices of the 
parameters. The stability and second-order convergence 
are demonstrated in Fig. |l|, which displays $ and v as 
functions of r along an outgoing null ray deep inside the 
black hole for various values of the grid parameter N. 

Figure ^ displays equi-spaced (in u) outgoing null rays 
(with constant values of u) in the rv-plane. The strong 
nonlinear dynamics is demonstrated by the rapid increase 
in the apparent horizon. All the rays which are not out- 
side the event horizon and escape to infinity, either ter- 
minate at r = within a finite lapse of advanced time v 
(type I) , or approach a finite limiting value of r as v — > oo 
of r at large values of v (type II) . 

The null portion of the singularity: Figure |^(A) shows 
the behavior of r tV along a type-II outgoing ray. At 
late times r „ oc v~ 2 . Figure ||(B) shows the behav- 
ior of $..„ along the same outgoing ray. At late times 
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FIG. 1: Behavior of $ and r along an outgoing ray (at fixed 
u) for various values of N. Upper panel (A): v as a function 
of r. Lower panel (B): $ as a function of r. In both panels 
dotted lines correspond to N = 10, dashed lines to N = 20, 
and solid lines to N = 40, and the data are taken for u — 21.9. 



oc v . This implies that along type-II rays r in- 
deed approaches a non-zero finite value as v — * oo, but 
4> diverges logarithmically in the same limit. This be- 
havior is in sharp contrast with the behavior of $ in the 
case of perturbations with compact support, where 
too, approaches a non-zero finite value. We next check 
the detailed behavior of the fields along type-II rays. We 
find that type-II rays terminate (in the infinite future as 
v — > oo) at a curvature singularity. This is demonstrated 
by Fig. §(C), which shows the exponential increase in 
R = (RapR 01 ^) 1 ^ 2 along the same outgoing ray, R a p be- 
ing the Ricci tensor. (Notice that here R is not the Ricci 
curvature scalar.) The finiteness of r at the singular hy- 
persurface suggests that the non-central portion of the 
singularity is deformationally weak. 

Note that the late-time behavior of r tV and as is 
clear from Fig. ^ starts to dominate much earlier than 
in the case of perturbations of compact support. In par- 
ticular, no quasi-normal modes (QNM) are visible. The 
reason for that is that the gradients of $ decays here 
much slower than the tails in the case of perturbations 
of compact support. Specifically, decays here accord- 
ing to an inverse power-law with a smaller index than 
in the case of perturbations with compact support. (In 
the latter case the index is 4 for spherical perturbations, 
whereas here we have an index of 1.) Because the field 
is stronger, it starts dominating earlier, and overwhelms 
the rapidly-decaying oscillations of the QNM. 

The spacelike portion of the singularity: Type-I rays 
terminate at a spacelike singularity. The spacelike sin- 
gularity inside a spherical charged black hole which is 
perturbed by a scalar field was studied within a simpli- 
fied homogeneous model in Ref. |l5j ] , where the pointwise 
behavior of the geometry and the field was found. It was 
also shown in p5| that approaching the spacetime sin- 
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FIG. 2: Outgoing null rays (with fixed values of u) in the rv- 
plane. The solid lines correspond to different fixed values of 
u, and the dotted line describes the apparent horizon, which 
approaches u = 21.07185 at late values of advanced time v. 
The values of u for which rays are shown (from right to left) 
are from u = 20.8 to 22.7 in equal increments of Ah = 0.1. 
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FIG. 3: Behavior of the fields and the curvature of spacetime 
along a type-II outgoing null ray. Upper panel (A): r >v as a 
function of v. Middle panel (B): $,„ as a function of v. Lower 
panel (C): R as a function of v. The data are shown along 
u = 21.1. 



gularity, the fully nonlinear and inhomogeneous numeri- 
cal solution (where the perturbation field had a compact 
support on the characteristic hypersurface) was in full 
agreement with the pointwise behavior. The study of the 
singularity in Rcf. |l5[ ] was local: no assumptions were 
made regarding the form of the perturbation field on the 
characteristic hypersurface. We thus expect that type-I 
rays terminate at a spacelike singularity whose pointwise 
behavior is well described by the singularity of Ref . p5| . 
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FIG. 4: The scalar field "I" and the radial coordinate r ap- 
proaching the singularity along a type-I null ray. Upper panel 
(A): d, - u as a function of r. Lower panel (B): $ as a func- 
tion of r. The data are shown along the outgoing null ray at 
u = 21.9. 

Assuming homogeneity, one finds approaching the space- 
like singularity inside a spherical charged black hole with 
a scalar field, that 

$(r) = v //3Tlliir + O(r^). (8) 

Here, (3 > is a constant (which numerically can be 
found to depend on u). Also, along an outgoing null ray 
one can show that, to the leading order in j;» - d, 

r(u Q ,v) = {d 2 /(f3+l)}V\v*-v) 1 / 2 , (9) 

where d is a gauge-dependent quatity (which depends on 
the scaling of the temporal coordinate), and = v(r = 
0) along that null ray. 

The agreement of our results with Eqs. (||) and @ is 
already apparent from Fig. |l|. Next, we check this agree- 
ment in greater detail. In Fig. ||(B) we show the behavior 
of the scalar field $ as a function of r(v) along a type-I 
outgoing null ray. This logarithmic behavior is consis- 
tent with Eq. (JsJ) . Along all type-I rays we find the same 
logarithmic divergence of <E>, including along rays which 
initially are outside or inside the apparent horizon. (The 
only difference between different rays is that the slope 
of the graph, i.e., the value of the parameter (3, changes 
from one ray to another.) Figure ^(A) displays — v 
as a function of r along the same outgoing null ray. The 
asymptotic behavior approaching the singularity agrees 
very nicely with Eq. (|9j). 

We conclude that the pointwise behavior at the singu- 
larity which we find in our simulations is well described 
by the singularity described in Ref. |l5| . Notice that this 
singularity is different from the Schwarzschild singular- 
ity: The former has (3 > 0, and the latter has (3 = — 1. 
This portion of the singularity then is scalar curvature, 
spacelike, and deformationally strong. 
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Conclusions: We studied the evolution of spacetime, 
and specifically the formation of curvature singularities, 
for a very simplied toy model of a spherical charged black 
hole, which is perturbed nonlincarly by a self-gravitating, 
spherical scalar field, which has non-compact support on 
the characteristic initial hypersurface. Although these 
perturbations are stronger than those which result from 
an initial profile with compact support (the gradient of 
the scalar field $ decays at late times as w _1 in our case, 
and as v~ 4 in the case of perturbations with compact sup- 
port) and consequently the evolution of spacetime and in 
particular of spacetime curvature is indeed dominated by 
the non-compact perturbations rather than by the per- 
turbations due to the collapse, we find that a portion of 
the Cauchy horizon still survives as a non-central, null 
singularity, rather than being utterly destroyed and re- 
placed by a central, spacelike singularity. The null gener- 
ators of the Cauchy horizon contract with retarded time 
u, and eventually arrive at r = 0, where the causal struc- 
ture and the strength the singularity change: the central 
singularity is spacelike and deformationally strong. This 
situation and the global causal structure is therefore very 
similar to that of a black hole perturbed by a pertur- 
bation field with compact support, despite the different 
details of the dynamics. 

The reason why the Cauchy horizon survived the in- 
troduction of non-compact perturbations as a null, non- 
central singularity is the following: We chose the charac- 
teristic field to be such, that although it is non-compact 
and does not belong to the same class of behavior on 
the event horizon at late advanced time as the tails, its 
gradient does. Specifically, along the event horizon $ is 



logarithmic in advanced time. This certainly does not 
belong to the class of the tails, which decay as an inverse 
integral power of advanced time on the event horizon. 
However, </> : „ decays like z; -1 along the event horizon at 
late advanced time, such that it does belong to the same 
class as the gradients of the tails. (Note, that no tails 
would ever produce n = 1, as n is at least 3 for all tails.) 
It is, in fact, <f>. v which is the important quantity, as cur- 
vature depends on the gradients of $, rather than on <f> 
itself. This particular form of jtI implies that 1Z can still 
be twice integrable approaching the Cauchy horizon. 

We therefore conclude that by themselves, perturba- 
tions of non-compact support, even when they dominate 
the dynamics, are not sufficient to obliterate the null, 
non-central singularity at the Cauchy horizon. It remains 
an open question, however, whether other classes of non- 
compact perturbations behave similarly. The CBR per- 
turbations, which are a generic source of perturbations 
for realistic black holes, are of particular interest. The 
energy influx of the CBR decays only on very long time 
scales due to the expansion of the universe (in a matter- 
dominated universe). (In a dark-energy dominated uni- 
verse the influx of CBR energy decays faster.) It is inter- 
esting to investigate how the Cauchy horizon singularity 
is affected by such perturbation fields, and also to inves- 
tigate which families of perturbing fields may destoy the 
null, non-central singularity. 
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